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Abstract 

The separation of a heavy quark and antiquark pair leads to the formation 
of a tube of flux, or string, which should break in the presence of light quark- 
antiquark pairs. This expected zero temperature phenomenon has proven 
elusive in simulations of lattice QCD. We present simulation results that show 

that the string does break in the confining phase at nonzero temperature. 
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In the absence of light quarks the heavy quark-antiquark potential is known quite accu¬ 
rately from numerical simulations of lattice quantum chromodynamics At large separa¬ 
tion i?, the potential rises linearly, as expected in a conhning theory. The potential in the 
presence of light quarks is less well determined because of the substantially higher computa¬ 
tional expense. Still, all the existing lattice data at zero temperature agree in that they 
do not show any indication of string breaking which would be signalled by a tendency of the 
potential to level off at large distances. The distances covered so far extend up to R<2 fm 
while it has been proposed that the dissociation threshold would be reached at separations 
somewhere between 1.5 and 1.8 fm P,^]. 

In this communication we present simulation results which show string breaking. These 
results have been obtained at nonzero temperature in the conhning phase of QCD. Our work 
conhrms trends found from simulations at nonzero temperature on smaller lattices 0. We 
have simulated QCD with two light havours of staggered dynamical quarks on lattices of size 
16^ X 4 (new work) and 12^ x 6 (conhgurations from Ref. ) at hxed values for the quark mass 
of mq/T = 0.15 and 0.075 respectively. The couplings were chosen to cover temperatures 
T below the critical temperature Tf. in the range of approximately O.TT^ < T < T^. The 
(temperature-dependent) heavy quark potential V{R,T) was extracted from Polyakov loop 
correlations 

(L(O)Lt(R)) = c exp{-R(|R|,T)/T} (1) 

where 

1 N.,-1 

L(f) =-tr n Uo(x,r) (2) 

T=0 

denotes the Polyakov loop at spatial coordinates x. In the limit R oo the correlation 
function should approach the cluster value |(L(0))p which vanishes if the potential is rising 
at large distances (conhnement) and which acquires a small nonzero value if the string 
breaks. 

In Figures ^ and ^ we present our data for the potential, at the values of (3 analyzed, in 
lattice units. The critical couplings (3c have been determined as 5.306 for W = 4 and 5.415 
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for Nr = 6 respectively. The Polyakov loop correlations have been compnted not only for on- 
axis separations bnt also for a conple of off-axis distance vectors R. Rotational invariance is 
reasonably well recovered if one nses the lattice Conlomb behavionr to determine the qnark- 
antiqnark separation, |i?| = l/Giat(R)- The data in Fignres ^ and Q qnite clearly show a 
flattening of the potential at lattice distances of abont 3 to 4 lattice spacings, depending on 
(3. Moreover, the height of the potential at these distances is in nice agreement already with 
the inhnite distance, clnster valne, shown as the right-most data point in each of the plots. 

In order to obtain a rongh estimate of the corresponding temperatnres in nnits of the 
critical temperatnre we applied the following procednre: at the given [3 and rriqa valnes 
an interpolation formnla [0] was ntilized to estimate the vector meson mass myo in lattice 
nnits as well as the ratio of psendoscalar to vector meson mass, mps/my. By means of a 
phenomenological formula which interpolates between the (experimentally measured) p and 
K* mass as function of the ratio mps/my, a value for my in physical units at the simulation 
quark mass can be obtained. This number is then used to estimate the value for the lattice 
spacing. This is certainly a rather crude procedure, yet, the resulting values for T/Tc = ad a 
show quite stable behaviour under variations of the procedure. The resulting temperature 
ratios are summarized in table |. 

Finally, in order to facilitate a comparison of the iV,- = 4 and 6 results with each other 
and with quenched data, the absolute scale was determined from a conventional Wilson loop 
measurement of the string tension at zero temperature at the critical /3c values. The Wilson 
loops did not show string breaking at the separations which could be explored. The results 
for the critical temperature in units of the string tension are obtained as Tdda = 0.436(8) 
for W = 4 and Td = 0.462(9) for W = 6 |§]. There is a substantial (25%) difference 
between scales set by my and a/^, which suggests a magnitude for systematic errors in the 
scale estimate. 

In Figure ^we show the potential in the presence of dynamical quarks in physical units. 
The data has been normalized to the cluster value i.e. the self energies have been taken 
out. The potential is flat within the error bars at distances larger than about 1 fm. It also 
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seems that the turn-over point is slightly T dependent, becoming smaller with increasing 
temperature. It is beyond the scope of the quality of the data at this stage, however, to 
quantify this statement. 

Assuming that the Wilson loop string tension is not affected by the absence of dynamical 
fermions one can then immediately compare quenched and full QCD potentials in physical 
units at the same temperature, as is shown in Fig. The quenched data has been taken 
from P] and was obtained in the same way, i.e. computed from Polyakov loop correlations. 
Each data set has slightly been shifted up or down to give rough agreement at interme¬ 
diate distances around 0.3 fm. Figure ^ contains, for further comparison, the dashed line 
denoting —7r/(12i?) -|- (420MeV)^i? which gives a good description of the zero temperature 
quenched potential. Note that the nonzero-temperature quenched potential is rising with 
distance R but the slope decreases with temperature, i.e. the (quenched) string tension is 
temperature dependent and becomes smaller closer to the critical Tc. Again, the comparison 
with quenched potentials at the same temperature demonstrates quite nicely that the po¬ 
tential in the presence of dynamical quarks becomes flat within the error bars at distances of 
about 1 fm. From Figure ^we conclude that the observed string breaking, albeit at nonzero 
temperature, is an effect caused by the presence of dynamical fermions. 

We have seen that string breaking is relatively easy to observe in the Polyakov loop 
correlation, while it is difficult to detect through the conventional Wilson loop observable. 
Why is this so? The Wilson loop observable creates a static quark-antiquark pair together 
with a flux tube joining them. In the presence of such a static pair at large R, we expect 
the correct ground state of the Hamiltonian to consist of two isolated heavy-light mesons, 
however. Such a state with an extra light dynamical quark pair has poor overlap with the 
flux-tube state, so it is presumably revealed only after evolution to a very large T. An 
improved Wilson-loop-style determination of the heavy quark potential in full QCD would 
employ a variational superposition of the flux-tube and two-heavy-meson states I30- 
The Polyakov loop approach, on the other hand, although limited in practical application 
to temperatures close to or above Tc, builds in no prejudices about the structure of the 
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static-pair ground state wave function. Screening from light quarks in the thermal ensemble 
occurs readily. 
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TABLES 


TABLE 1. Estimates of the temperature at the various (5 and Nj- values. 


II 

/3 

5.10 

5.20 

5.25 

5.28 


TIT, 

0.67 

0.79 

0.87 

0.94 

Nr = 6 

a 

5.37 

5.38 

5.39 

5.41 


TjT, 

0.85 

0.88 

0.91 

0.98 


7 





















0 0.5 1 1.5 2 2.5 3 3.5 4 


FIG. 3. The potential in physical units at various temperatures. The results are from lattices 
with Nt = 4 and 6, indicated by the number in brackets. The data has been normalized to the 
cluster value. 
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FIG. 4. The potential in physical units at various temperatures. Compared are quenched (open 

symbols) and full (filled symbols) QCD potentials at the same temperature. The dashed line is the 
zero temperature quenched potential. The data has been normalized to agree at distances around 
0.3 fm. 
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